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Abstract 

We consider super-inflating solutions in modified gravity for several popular families 
of f{R) functions. Using scalar field reformulation of /(i?)-gravity we describe how the 
form of effective scalar field potential can be used for explaining existence of stable 
super-inflation solutions in the theory under consideration. Several new solutions of 
this type have been found analytically and checked numerically. 

1 Introduction 

Theories of modified gravity first introduced in order to match General relativity and quan- 
tum physics, became very popular recently after discovery of accelerated expansion of our 
Universe. Being higher order, equation of motion of these theories can contain some regimes 
impossible in General relativity with its 2-d order equations. Some of new regimes may be 
useful for explaining observed accelerated expansion, others can spoil standard cosmology, 
so constructing viable models appears to be rather tricky (however, not impossible). 

On the other hand, though modification of an early Universe dynamics gives one of the 
earhest proposed inflation scenario (so called Starobinsky inflation, [1]), it is less popular now 
due to success of scalar field inflation theory [2], [3]. This theory does not require any modifi- 
cation of General Relativity, so any correction to high-energy dynamics being pre-infiational 
have no consequences today. However, this statement is valid only if such corrections do 
not modify significantly inflation itself. In particular, stable super-inflating vacuum solutions 
(solutions with increasing Hubble parameter) are dangerous. Recent analysis shows that sta- 
ble phantom solutions are common in f{R) gravity with power-law f{R) [1],[5]. This regime 
leads to Big Rip singularity in a finite time interval [6] [7], and can not be followed by some 
other non-singular regime. 
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Apart from the phantom, other super-inflating solution in i? + aB? gravity in known 
for decades. It have been discovered initially in [8j and rediscovered later ([9j[10]). The 
main feature of this solution is that the Hubble parameter grows linearly in time, so it takes 
infinite time before it reaches infinity. In generalized modified gravity other forms of this 
solution are possible, in particular RUR corrections in the action lead to ~ t^^'^ solution 
[TT] . Moreover, the presence of matter with some exotic equations of state or rather exotic 
functions f{R) leads to super-inflation with H [12] , [I3j , , ^15] , [16] in these papers such 
solution is called as " Little Rip" . However, from physical point of view, Ruzmaikin solution 
with its analogs and phantom solutions are not substantially different - both lead to Planck 
density reached in finite time interval, and prolongation of the solution further in time has 
no any physical justification. Both solution, if realised in early epoch of Universe evolution 
can precede inflation only being unstable. It should be noted that recently some models 
describing "graceful exit" from a phantom regime have been constructed [TH |T7], however, 
they are dealing with a phantom matter as a source of this regime, and can not be directly 
applied to /(-R)-cosmology. 

In the present paper we try to explain why stable super-inflating solution are so natural 
in f{R) gravity and present some of them which have not been found earlier. 



2 Effective scalar degree of freedom 



f{R) gravity is known to be dynamically equivalent to the General Relativity with the min- 
imally coupled scalar field ([H], [19], [20], [21]). Therefore, we expect that properties of cos- 
mological solutions in /(i?)-gravity can be guessed from the behavior of effective scalar field 
potential. We begin with the equations of motion and their trace [21] for /(i?)-gravity in the 
Jordan frame . 



f'G,.-f;^,+ [Df'-^{f-f'R) 



T, 



flU 1 

T 



Df' = -{2f-fR)+^ . 
Following [22], we introduce the scalar degree of freedom and define its potential: 



V'i 



= /'-2, 
^l(2/-/'i?) 



1 1 (2/ - f'R)^^dR =\j{V- f'R)f"dR. 



Then the trace equation becomes: 

""^The signature of the metric is assumed to be (— , +, +, +) and c = SttG = 1. 



(1) 
(2) 

(3) 
(4) 

(5) 
(6) 

(7) 
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Finally, in the fiat FLRW metric (rfs^ = —dt'^ + a(t)^dx^) the above equation is 

= -r'(0) + J^-3i/0. (8) 

After the substitution — > x we obtain a simple one- dimensional equation of motion for a 
"particle" in the potential V{x) with the friction ?)Hx and the driving force J-". This analogue 
is useful because behavior of the mechanical system is quite evident especially in the absence 
of external force. Using the Eqs.([5]),(|n]) one can see that the function 1^(0) is represented by 
a plane curve defined parametrically [0(-R), V{R)], where R = R(t). We consider the vacuum 
case (T = —3J^ = 0) - this case corresponds to the absence of driving force. 
The equations of motion for FLRW metric have the form 

3fH' = ^{f'R-f)-3Hf', (9) 
R = 6{2H^ + H). (10) 

In the present paper we consider mostly regular functions f{R) which remain finite for 
finite R with its first derivativ^. As we will see further, properties of the effective potential 
allow us to find some asymptotic cosmological solutions for particular functions f{R), which 
appears to be useful for detection of a possible super-inflating behavior. 

Movement in the plane [4>,V] corresponds to the solution R{t). On the other hand, the 
composite function obeys the equation of motion for "particle" in classical mechanics 

([8]). Thus let us consider the "particle" motion. The following situations are possible: 

1. V{(j)) has a local minimum, hence in the limit of late time the "particle" would fall down 
in this minimum due to the Hubble friction. Scalar degree of freedom takes the value 
= (pa, for which V'icpa) = 0. Therefore, 2/(/?J - f{Ra)Ra = (see (HD) and we see 
that Ra = R{4>a) corresponds to the stable de Sitter (dS) {H = const), or Minkowski 
{H = 0) asymptotic solution. 

2. V{(f)) has neither local nor global minima, hence it is unbounded from below. The "par- 
ticle" must fall into an unrestricted area of the potential curve. If |^'(0)| is unbounded 
for — > — oo this, according to Eq.(j4]) indicates an infinite growth of |-R|. Stable super- 
inflating solutions in the signature chosen correspond to i? — )■ +oo solutions. The case 
of limited V{(j)) for infinite R should be excluded from our considerations. We have not 
seen it in particular families of f{R) functions studied in the present paper (see below). 

If V^(0) has the local minimum (possibly more than one) and is unbounded from below then 
different asymptotic regimes are possible depending on initial conditions. 

It should be noted that due to presence of non-standard singularities in f{R) gravity 
(which occurs when a coefficient at the highest derivative term in the equations of motion 
vanishes) the cosmological evolution does not necessarily ends up at some smooth asymptotic 
regime. That is why existence of such regimes can only be guessed from the form of effective 
potential and should be checked explicitly. In the present paper we write down analytic from 

^It should be also noted that for this theory to be perturbatively stable, one needs /' > to avoid ghosts 
([13) and /" > to avoid tachyon instability ([ ^ .[25 ] ). 



3 



of super-inflating solutions in question and check their stability only numerically, leaving full 
analytic treatment of stability (regarding solution which have not been known earlier) for a 
future work. As for possible unstable solutions (which may be interesting from the viewpoint 
of the Starobinsky inflation) the obvious example is a potential which has a local maximum. 
We will not consider unstable solutions further in this paper. 

It also should be noted that the proposed effective potential does not help us to study 
the regimes where f'{R) becomes zero at some finite R and the effective gravitational con- 
stant changes sign. To be sure that a phase curve can safely cross this point, one should 
apply the commonly used conformal transformation to the Einstein frame ( [18] , [19] , [20] , [21] ) . 
However, as the purpose of this work is to study super-inflating regimes in f{R) gravity the 
method proposed in [22] provides relatively simple scheme of detecting a super-inflation be- 
havior based on mechanical analogue because the dynamics of the effective scalar field in this 
parametrisation resembles particle motion with friction in the potential field. The standard 
way [19] may be useful in describing a global picture of corresponding dynamics which is 
beyond the scope of the present paper. 

3 Examples of /(i?)-functions and numerical results 

3.1 i{R)=R + aR^ 

In this case: 

V(0) = 3(2^"_^) (^'^(2iV' - 3iV + 1) - a{N -2){N- 1)NR'''-') , (11) 
= aNR^-^ - 1 . (12) 
Two examples shown in FigsIH 12 are 

V^(0) = ^(0' + 20) for f{R) = R + R\ (13) 
V{(t>) = -4 + 2(20 - 5)(0 + 1) for f{R)=R + R\ (14) 

00 



The case of N> 2. For > 2 (A^ = 4 on the figure, the local minimum is located at i? = 0) 
the potential has the asymptotics — )■ oo at i? — )■ — oo and V — )■ — oo at i? — )■ oo. There is 
no global minimum and the "particle" falls down into the region of large negative V. This 
leads to V — oo and, correspondingly, R — t- ocE|. The numerical procedure confirms our 
reasoning (see Fig|3]). Hence in the case N > 2 the curvature diverges and we obtain stable 

2N^ -3N+1 

super-inflating solution. This solution appears to be a power-law phantom a ~ t 2-iv and 
have been found earlier in [26], [27] [1]. Note, that the shape of the potential and asymptotic 
behavior does not depend essentially on the sign of a in this regime (see FigJl]). 

^Consider theory with odd N we have to choose one of the branches (i? > or i? < 0) because of the 
potential ambiguity. Also i? > for theories with non-integer TV. 
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Figure 1: Effective potential of 
scalaron for f{R) = R + R^. 



H 



Figure 3: Stable pliantom solu- 
tion for /(-R) = R + aR^ does 
not depend on tfie coupling a. 



Figure 2: Effective potential of 
scalaron for f{R) = R + R^. 




Figure 4: Effective potential of 
scalaron for f{R) = R — R'^. 



The case of N= 2. Effective potential of the scalar degree of freedom has the global 
minimum at R = for a > (see Figl2]) and the potential is unbounded from below and has 
asymptotic V — oo as R ^ ±oo for a < (see Figj5]). Therefore super- inflation (presented 
here by stable Ruzmaikin solution) exists in the last case (see Figj6]) and the case of positive 
a provides Minkowski solution (see Figs|7]). The stability of Minkowski solution might be 
seen only numerically because /'(O) ^ 0, /"(O) = 0, /'" = in this case (see [28]). Comparing 
with the previous case we note that changing the sign of the coupling constant in modified 
gravity with R + aR^ makes the super-inflation behavior stable. On the contrary, it does not 
give new result for R + aR'^ , N > 2. 



The case of 1 <N< 2. By the definition, the curvature in such case must be positive. The 
effective potential is bounded from below for any couplings (see FigJH]). The numerical study 
shows that the theory is free from super-inflating solutions in such case. Detailed analysis 
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of cosmological dynamics in this case was done in works [H [5] and it is out of topic of our 
work. Also we do not consider the case < 1 through the paper because it is sub-dominant 
in the high-curvature regime and, therefore, not essential for super-inflation. 




Figure 5: Effective potential of Figure 6: Stable super-inflating 

scalaron for f{R) = R — R^. behaviour for f{R) = R — R^. 




The results of the numerical procedure are in good agreement with those obtained by using 
the Dynamical system approach in [4j,[5]. In fact, we have only phantom stable power-law 
solutions for N > 2. 
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3.2 f(R)=i?^exp(i?) 
In this case: 

2R 

y{<t>) = -^[N{N -1){N- 2)(T2jv-2 + (3iV2 - bN)(T2N-i + (3Ar - 2)a2N + , (15) 



V- (-1)' k\ 



1=0 



(16) 
(17) 





Figure 9: Effective potential 
of scalar degree of freedom for 
f{R) = R^exp{R), 1<N <2. 

\ 



H 



Figure 11: Stable asymptotic so- 
lution for f{R) = Re^ corre- 
sponding to the local minimum 
point at effective potential. 



Figure 10: Effective potential 
of scalar degree of freedom for 
f{R) = i?^exp(i?), N>2. 



H 




Figure 12: Stable asymptotic 
super-inflation solution for 
f{R)^ R^e^,N >1. 
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The potential has the asymptotic V — )■ —oo for </>—)■ oo giving R — )■ oo, hence the 
"particle" falls down in the potential hole (see FigsJHl fTO]) - Moreover, the potential for 
f{R) = exp {R), 1 < N < 2 has the local minimum aX R = and depending on initial 
conditions, one can obtain the stable Minkowski or super-inflating solutions. The numerical 
procedure confirms our reasoning (see FigsJTT l [T^ . 

Let us explain obtained results analytically. At first, consider de Sitter solutions in the 
theory: 

/'(i?o)i?o-2/(i?o) =0,^ Ro = or Ro = 2 - N . 

Note that i? > for dS in the chosen signature. Stability condition for de Sitter solution [2B] 
requires < f"{Ro)Ro < f'{Ro): 

fo,^.0: ».^^<1 fo.. N<2. (19) 

J [Ro) 

Indeed, Minkowski solution is stable for 1 < < 2 and de Sitter solution is not. 

Secondly, diverging numerical solution demonstrates linear behavior of the Hubble pa- 
rameter hence one could expect that Ruzmaikin-like solution exists in such theory. Consider: 

a{t) = aoe fe+i , ^^g) 
H{t) = Bt'' . 

Substituting this ansatz in the equation of motion ([9]) for f{R) = R^ exp (R), we find: 

+ Ast"^' + A^e^+^ + Aiot''^ + Ant''^+' = . ^'^^^ 

The biggest power indexes (which we should equate) belong to the terms A^t^^^^ and Ant^^^'^. 
This means that 

4A; + 4 = 5A; + 3 , 
k = l. 

Now we find the constant B. One should set the sum of coefficients in front of dominating 
term A^ + An to zero: 

576A;5^6^ - 2AB%^ = , ^ = 1/24 . (22) 

Thus the Ruzmaikin-like solution a ~ exp (t^/48) exists in the modified gravity theory with 
exponential correction for any > 1. Note, that the Einstein term in the case of f{R) = 
R + aR^e^ is sub-dominant, and, therefore, is not necessary for Ruzmaikin solution here in 
contrast to the R + R^ case (there is no such solution in a pure R^ gravity). After neglecting 
the Einstein term the constant a drops down from the equations of motion, and, therefore, 
properties of the solution does not depend on the sign of a. 

This theory have been studied using dynamical system approach in |4] , though Ruzmaikin 
solution have not been found there, so the list of asymptotic solutions from |3] is incomplete. 
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3.3 f{R)=R + aR^lnR 

This function has some interesting features. We could expect that the logarithmic correction 
is not essential for high-curvature regime, and this theory would be equivalent to f{R) = 
R + aR^ described above in this paper. However,this correction allow us to obtain new 
solutions for = 1, 2 which do not exist without them. It is worth to note that despite 
logarithm diverges at zero, the chosen form of f{R) ensures that it is regular in this point. 
What can happen for a singular f{R), we discuss in the next section. 

Let us write down the equations for scalar field and its potential in the terms of the R as 
usual. 



\3 



(-16 N'^ -UN^ + 2N + 4:N^ + 21 N^) aR^^'^ (In Rf 

(23) 



3(2iV- 1 

+ (21 _ ^0 iV + 2 + 8 A^^ - 22 N'^)aR^^-^ In R 
+ {20 - 1 - 8 - 18 N"^ + 7 N)aR^ \nR 



+ {AN^ -3N'^ + l)aR^^-^ + {1 + 12 N"^ - 8 - Q N)R^ 

(t)=a{N\nR+l)R^-^ -1. (24) 
For example in the case of = 1, a = 1 these expressions reduce to: 

V{(p) = -R{]xiR-l), (j) = lnR. (25) 
3 

The case of N> 2. For all > 2 the effective potential has no global minimum. Thus 
this theory should contain regimes with diverging curvature. Numerical calculation shows 
stable phantom solution in the considered case (see FigsJTSl and [T^ . The form of phantom 
solutions found for f{R) = R + aR^ in [1] does not change (it is easy to see using power-law 
ansatz), which is natural to expect due to very slow growth of the logarithm function: 

X 2JV^-3iV+l 

a{t) = ao{t - to) . (26) 

Note that the numerical solution does not depend on the coupling in front of logarithmic 
term as well as the shape of effective potential in this case. 



The case of N= 2. We remind a reader that for f{R) = R + aR^ the case of A^ = 2 
only gives us the Ruzmaikin solution (even if we allow non-integer power indexes). Effective 
potentials in the theory with logarithmic correction for N = 2 and A^ > 2 do not depend 
much on coupling a and also do not significantly differ from each other (see FigsJT3| [T5] and 
[TBI) . All such potentials are unbounded from below, however, the Ruzmaikin solution existing 
in the theory without logarithmic correction, disappears when we take this corrections into 
account. Choosing the ansatz in the form H = Bt'^, we obtain for f{R) = R + aR^ InR: 

f + Ait''+^ In {UBH^'' + QBkt^-^) + Ait^^^'^ + A3 In {UBH^'' + 6Bkt''~^) + A^t'^^^ + A5 = . 

(27) 
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Figure 13: Effective potential for 
f{R) = R + R^\nR, N>2. 



H 



Figure 14: Stable phantom so- 
lution for f{R) = R + R^lnR, 
N >2. 



The first term originates from the Einstein term. The third term becomes dominant for large 
times and it can not be compensated by another terms. Therefore, there is no Ruzmaikin 
solution in this theory. Using the power-law ansatz a{t) = aQ{t — to)^ one can easily obtain 
that there is no simple power-law phantom solution too. However, our numerical procedure 
shows the super-inflating solution diverging at a finite time as presented in Fig. [171 

Let us explain the nature of such solution in the theory with logarithmic correction. 
Consider: 

Hit) = a'"l^('°-'>1 = + , (28) 

to — t to — t 

Substituting this ansatz in the equation of motion ([9]) we find: 

- 365 \n{R) ln^[A(to - t)] + 6 = . 

(29) 



^(to - tf ln2[A(to - t)] + 42 ln[A(to - t)] + (545 - 30) \n^[^{tQ - t)] - hAB\n^[A{tQ - t)] 
5^ ln^[A(to - i)] + 6 ln(i?) + 24 ln(i?) ln[A(to - t)] + (365 - 18) ln(i?) ln^[A(to - t)] 



In the high-curvature case 

1252 In^ [A(to - t)] 65 65 ln[A(to - t)] 

ln(i?) ^ -21n(to -t) + 0(lnln | to - ^ I) , (30) 

so terms with 125^ and —365 factors dominate {t — )■ to) and lead to approximate solution 
H ~ ln[A(to — t)]/(to — t), which have been found numerically. The A constant determines 
the time scale for which this solution exists and depends on the initial conditions. To obtain 
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Figure 15: Effective poten- 
tial for f{R) = R- R^lnR, 
N>2. 



Figure 16: Effective potential 
for f{R)^R + R^lnR. 



H 




Figure 17: Stable super- 
inffating solution for f{R) — 
R + R^ InR. 



Figure 18: Effective potential 
for f{R) ^R + RlnR. 



H 





Figure 19: Stable de Sitter so- 
lution for f{R) ^R + RlnR. 



Figure 20: Effective potential 
for f{R) ^ R- RhiR. 
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the constant B one should set the coefficient in front of dominating terms to zero (factor —2 
arises from f l30|) ): 

125^ + 725 = ^ B = -6 . (31) 

The Einstein term is not important for this regime in such theory as well as in the previous 
subsection, thus stable super-inflating solution exists for any a 7^ 0. 

The case of 1 < N< 2 The potential has a global minimum (at = 0, i? = 1 in the 
case a = 1,A^ = 1), fiat small curvature asymptotic — )■ —00 corresponding to i? — )■ 
and exponential potential wall for large curvature (see Fig JTH]) . Thus this theory does not 
contain solutions with diverging curvature. The numerical procedure demonstrates stable 
de Sitter regime (see Fig JT9|) in the case 1 < < 2 for each certain coupling a even for N 
close to 2. Effective potentials in this case does not depend significantly on the sign in front 
of logarithmic term (see FigsJTKl and [2U]) . The existence and stability of de Sitter solution 
(as well as the absence of super-inflating solutions) for = 1 have been proved using the 
dynamical systems methods in [22] . 

Let us study the existence and stability of de Sitter solution in the case iV = 2 — e,e^ 1: 

nR)R-2f{R) = 0,^ (32) 
R + a{N \n{R) + 1)R^ -2R- 2aR^ \n{R) = , (33) 
-R + aR^[{N\n{R) + 1) -21n(i?)] =0, (34) 

The numerical procedure shows existence of a stable de Sitter solution with large curvature 
-R ^ 1. Thus the Einstein term is negligible for e ^ 1 and a drops out from the equations 
of motion: 

aR^[{N -2)\n{R) + l]=0, (35) 
R = ei/(2-^) . (36) 

Stability condition < f"{Ro)Ro < f'{Ro) (see [2B]) in the considered case implies: 

2 - 72 < < 2 . (37) 

Indeed, in the theory with R + aR^~^ In R we have no super-inflating solutions and obtain 
stable de Sitter solution with Rq f» e^/^. Note that both solutions ( (1261) and ( l36l) ) with 
A'^ = 2 ± e are sufficiently non-perturbative and limiting e — )■ does not allow us obtain 
solution existing only for N = 2. The situation here is the same with the known case 
R + ai?^"*"^, where Ruzmaikin solution exists only for e = and does not exist for any 
non-zero e. 

As in the former case, the de Sitter solution is stable, the only important difference is 
that the existence and stability of de Sitter solution with the logarithmic corrections does 
not depend upon the sign of the coupling constant a in the limit of small e. In the general 
case of an arbitrary A^ G [1, 2] the situation depends on both A^ and a, detailed description 
of de Sitter solution for this case is beyond the scope of this paper. 
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3.4 i{R)=R-l/R 

In this subsection we discuss that can happen if f{R) and (or) V{(l)) is singular for the finite 
value of its argument. In this case unboundedness of the effective potential from below does 
not indicate a super-inflating solution because R can tend to a finite limiting value for which 
f{R) diverges. Consider the prototype of General Relativity infra-red modification, proposed 
in [29] . It has been already studied that the scalar degree of freedom is unstable in this theory 
(|25]). but it is illustrative example in which mechanical analogue still works and gives results 
different from obtained above. The scalar field for this model is: 



2 



3 i?3 

1 



(38) 
(39) 



V{(f)) is an ambiguous function of 0, thus consider the branch R> (see Figs J2T] and l22l). 




H 




Figure 21: Effective potential for 
f{R) = R-1/R, R>0. 



Figure 22: Stable numerical so- 
lution for f{R) = R- 1/R. 



The potential has the asymptotic V — )• — oo for — )■ oo corresponding to i? — )■ , so the 
"particle" falls down in the potential hole and provides a stable Minkowski solutioiQ. As a 
result, the potential unbounded from below does not lead to a super-inflating solution in this 



case. 



The stability was checked only numerically because = (see [5F) in this case. 
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4 Conclusion 



It this paper we have explored the idea of determining asymptotic cosmological behavior in 
/(i?)-gravity using the effective potential of scalar degree of freedom for regular functions 
f{R). The presence of the global minimum appears to be incompatible with super- inflating 
cosmological solutions, on the other hand, unboundedncss of the potential from below indi- 
cates existence of such solutions. We considered only vacuum case; in the presence of matter 
the dynamics is not so evident. 

As an application of this property we considered three families of functions f{R) very 
popular in current researches on modified gravity, namely R + aR^ , R + aR^ exp{R) and 
R + aR^ In R. It appears that only the following potentials from these infinite families are 
free from stable super- inflation solution: the known R + aR^ , 1 < N < 2; R + aR^, a > 
cases and the case oi R + aR^ InR ,1 < N < 2. All other f{R) functions from these famihes 
contain stable super-inflating solutions, some of them have not been discovered earlier. 
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